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ABSTRACT 


This  project  was  a  theoretical  investigation  of  particulate 
impaction  on  spheres  and  cylinders.  The  motion  model  developed  was 
Implemented  on  a  computer  and  yielded  results  focused  on  two  main  goals: 
first,  the  net  effect  of  gravity  on  particulate  impaction  was  determined; 
and  second,  a  man  simulation  was  conducted.  This  simulation  calculated  to 
a  first  approximation  the  amount  of  chemical  that  would  impact  on  a  man 
subjected  to  a  chemical  attack. , 
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INTRODUCTION 


1.  In  assessing  the  hazard  to  troops  from  attacks  with  chemical 

agents  in  aerosol  form,  it  Is  important  to  be  able  to  predict  the  motion  of 
the  chemical  agents  as  they  disperse  through  the  atmosphere  and  impact  on 
various  target  objects  (Figure  1).  This  project  was  conducted  to 
investigate  the  impaction  aspect  of  the  overall  motion  problem  specifically 
to  provide  an  answer  to  the  question,  "How  much  chemical  will  hit  a 
target. 
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2.  Some  theoretical  studies  have  been  performed  in  this  area  of 
aerosol  impaction1.  This  project  sought  to  upgrade  the  theoretical  models 
used  through  the  inclusion  of  the  gravity  force  into  the  problem;  in  fact, 
the  secondary  purpose  of  the  project  was  to  quantify  the  error  incurred 
when  gravitational  effects  are  ignored.  The  primary  purpose  of  the  project 
was  to  develop  a  model  which  would  calculate  the  amount  of  chemical  that  a 
man  received  in  relation  to  how  much  the  ground  received. 

3.  The  complex  nature  of  this  fluid  mechanical  problem  required 

a  number  of  simplifying  approximations  in  order  to  be  tractable. 
Consequently,  the  results  obtained  are  first  approximations  only, 
calculated  on  the  basis  of  a  theoretical  model  which  required  computer 
programs  for  implementation.  The  project  results  are  seen  as  a 
foundation  upon  which  to  conduct  experimental  and/or  improved  theoretical 
research. 


4.  The  project  evolved  significantly  durings  its  lifetime.  The 

initial  phases  were  very  much  a  learning  phase  during  which  the  model  was 
continuously  tested  and  changed.  Spherical  targets  were  used  for  these 
initial  tests,  tests  which  also  provided  much  data  on  the  effect  of  gravity 
on  the  problem.  This  initial  phase  ended  with  the  study  of  cylindrical 
targets  during  which  the  last  gravity  effect  calculations  were  performed. 
The  final  phase  of  the  project  dealt  with  the  man-target  simulation. 

BACKGROUND  THEORY 

Fundamentals  of  Particulate  Motion 


5.  The  derivation  of  the  aerosol  particle  motion  equation  follows 

directly  from  Newton's  Second  Law2.  There  are  two  forces  present,  namely 
gravity  and  aerodynamic  drag: 
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*  >1)  -  v|2>  [1] 

Ii="2  (2] 

The  particles  are  assumed  to  be  spherical  because  of  their  small  size 
(<2000  un  radius).  The  resulting  equation  of  motion  Is  therefore: 

mv  =  C0(itrp2)(ispa  |u  -  V  j2)  +  m£  [3] 


Noting  that; 


m 


[4] 


Re  = 


2rppa 


u  .  V 


vie  can  re-arrange  equation  [3]  to  produce: 


3uC0Re(U  -  V) 


16r_2p_ 

P  P 


[5] 


[6] 


6.  Equation  [6]  can  be  non-dimensional  1  zed  by  using  the  free  stream 

velocity  U  and  the  characteristic  target  length  L: 


v.'  =  U_lv 

[7] 

•>l 

CM 

1 

=> 

—1 

II 

•>| 

[81 

u'  =  U" *0 

[9] 

t'  =  tUL"1 

[lOj 

g'  -  LU-2£ 

HI] 
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If  we  define: 


K  = 


2p„r „2U 
p  p 


9uL 


[12] 


U2 

F  =  [13 J 

Lg 

then  equation  [6]  can  be  written  in  the  following  non-dimensional  form: 


CDRe{U'-v')  1 

-  +  _ 

24  K  F 


[14] 


7.  Clearly,  the  problem  can  be  solved  by  either  equation  [6]  or 

equation  [14].  Traditionally,  these  problems  have  been  solved  in  non- 
dimensional  form  resulting  in  graphs  of  collection  efficiency  vs.  inertia 
parameter  for  various  Reynold's  numbers  (Collection  efficiency  is  explained 
in  the  next  section).  In  this  problem,  however,  the  inclusion  of  gravity 
adds  another  parameter,  namely  the  Froude  number.  The  graphical  presenta¬ 
tion  of  non-dimensional  results  now  becomes  quite  complicated,  and  the 
physical  interpretation  of  such  results  becomes  obscure.  It  was  felt  that 
superior  physical  insight  and  applicability  would  result  from  analysis  of 
the  restricted  case  of  motion  in  air  under  representative  experimental 
situations.  Equation  [6]  was  therefore  used  in  the  model. 


8.  In  reference  frame  1  (Figure  2),  equation  [6]  can  be  resolved 

into  the  following  scalar  equations: 


\  -  bCDRe  (Ux  -  V 

[15] 

vy  =  bC0Re  (Uy  -  vy) 

[16] 

=  bCQRe  (Uz  -  vz)  -  g 

[17] 
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where 

b  = 


3w 


16r„*-p 


2  , 

P  P 


[18] 


9.  The  drag  coefficient  for  spheres  Is  defined  In  terms  of  the 

Reynold's  number  by  the  following  equations3: 

For  Re<4,  CgRe2 

Re  =  -  2. 3363x10" 4 (CnRe2)+2. 0154x10“ 6(CnRe2) 3 

24  u  u 

-  6.9105xl0"9(CDRe2)4  [19] 


For  3<Re<104,  log1Q  Re  =  1.29536  +  9.86  x  10"Mlog10CDRe2)-  4.6677 

x  10“2  ( 1 og  1Q  CpRe2)  +  1.1235  x  10"3  (log10  CpRe2)3  [20 ] 

10.  Two  assumptions  were  made  regarding  the  initial  condition  of  the 

aerosol  particle.  First,  It  was  assumed  that  the  particle  was  moving 
horizontally  at  the  free  stream  velocity  and  that  a  starting  position 
upstream  would  experience  negligible  flow  perturbations  caused  by  the 
target.  Second,  the  particle  was  defined  to  be  falling  at  Its  terminal 
velocity,  because  small  aerosol  particles  quickly  attain  that  speed. 
Consequently,  the  particle  Initially  possessed  a  velocity  given  by  its 
terminal  velocity  (v^  )  and  the  free  stream  velocity,  with  a  direction 
defined  by  the  angle  y  £o  the  x-axIs: 


y  = 


[21] 


The  terminal  velocity  was  calculated  by  setting  v  =  0  and  Uz  =  0  in 

equation  [17],  and  solving  for  v7  .  The  result  was: 

^  o 
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3paCD 


[22] 


11.  Therefore,  simultaneous  solution  of  [5],  [22]  and  one  of  [19]  or 

[20]  yielded  the  terminal  velocity. 


12.  The  flow  equations  will  be  derived  In  paragraphs  15-24  for  the 

various  target  geometries.  All  such  equations  were  inviscid  fluid 
equations.  These  were  deemed  applicable  to  the  problem  because  real  fluid 
flow  Is  virtually  Ideal  on  the  leading  face  of  an  object,  which  was  the 
only  face  of  Interest  In  this  particle  Impaction  problem.  It  was  also 
assumed  that  the  chemical  particles  themselves  did  not  perturb  the  fluid 
flow. 


Collection  Efficiency 

13.  The  concept  of  a  collection  efficiency  Is  the  primary  means  of 
quantifying  the  amount  of  chemical  which  Impacts  on  a  target  (Figure  3). 
The  collection  efficiency  Is  a  number  ranging  from  zero  to  one,  and  is 
defined  by  equation  [23]: 

Cross-sectional  area  of  the  envelope 

Collection  = - 

Efficiency  Cross-sectional  area  of  the  target 
The  envelope  Is  that  region  on  the  starting  plane  In  which  initial  particle 

placement  results  In  a  trajectory  which  hits  the  target.  Clearly,  the 

absence  of  the  fluid  medium  would  result  In  a  collection  efficiency  of  one, 

because  the  particles  would  travel  in  straight  lines,  and  massless 

particles  would  have  a  collection  efficiency  of  zero,  because  the  particles 

would  follow  fluid  streamlines  around  the  target  without  Impaction. 

14.  The  task  at  hand  is  to  calculate  the  collection  efficiency  for 
any  given  test  situation.  This  Is  done  using  a  half-interval  method  to 
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find  the  boundaries  of  the  envelope.  A  series  of  particles  are  considered, 
each  of  which  begins  on  the  starting  plane  with  the  same  initial  velocity 
(Figure  4).  Two  initial  particle  positions  are  required  to  start  the 

procedure:  an  'inside'  one  that  results  in  target  impaction  (P2)  and  an 

'outside'  one  that  results  in  target  miss  (PI).  A  third  position  (P3)  is 
calculated  such  that  it  lies  halfway  between  the  first  two  points;  this 
particle  is  then  tracked  to  the  target.  If  it  misses,  the  (P3)  becomes  the 
new  'outside'  position  and  a  new  starting  position  is  chosen  at  (P4).  If 
it  hits,  then  (P3)  becomes  the  new  'inside'  position,  and  a  new  starting 

position  is  chosen  at  (P5).  The  distance  between  the  inside  and  outside 
positions  is  successively  reduced  by  half,  and  it  asymptotes  to  the 

envelope  boundary.  In  this  manner,  the  envelope  boundary  is  defined, 
allowing  its  area  to  be  calculated  for  use  in  equation  [23].  Note  that  the 
boundary  is  generally  a  two-dimensional  curve  on  the  starting  plane  which 
requires  a  large  number  of  boundary  points  before  an  accurate  estimation  of 
its  shape  and  size  can  be  made. 


Target  Geometry  For  Spheres 


15.  The  fluid  velocity  field  around  a  sphere  is  given  in  terms  of 

polar  coordinates4: 


Ur  =  U  cos  0 


[24] 


UQ  ■  -U  sin  0 

o 


1  +  1/2 


[25] 


These  equations  are  valid  for  any  plane  passing  through  the  origin  of  the 
sphere  and  parallel  to  the  x-axis  of  Frame  1.  Let  us  designate  such  a 
plane  as  Frame  0  with  x,  n,  axes;  Frame  1  is  related  to  Frame  0  by  the 


% 
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rotation  angle  0  about  the  x-axIs  (Figure  5).  In  terms  of  Frame  0 
cartesian  coordinates,  equations  [24]  and  [25]  become: 


A3 

Ux  =  Ucos2  0 

Hu 

+  Usin2  0 

‘ * u%)_ 

[26] 


-3U  /l\3 

U  =  -  cos  0  sin  0  [27] 

n  2  (0 


Equation  [28]  Is  obtained  by  Inspection.  These  three  equations  can  be  used 
to  specify  the  velocity  field  at  any  point  in  space  in  Frame  1  via  a 
suitable  rotation  of  Frame  0  to  align  It  with  the  position  of  Interest.  In 
matrix  notation,  the  required  velocity  components  are  generated  as 
follows: 


U„ 

1 

u* 

II 

< 

1  0  0 

0  cos0  -  sin3 
0  s1n0  cos0 


[29] 


Upon  evaluation  we  get: 


Ux  =  U  cos2© 


/L  3 

L  3 

-b  ] 

+  U  s1n2e 

‘ +  Vj_ 

[30] 


-3U/L  3 

U  =  I  cos  0  sin  0  cos  0 

y  2  lr 


[31] 
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-3U/L\ 3 

Uz  =  -  cos  O  sin  0  sin  e  [32] 

By  definition,  we  obtain  these  auxiliary  equations: 


x  -  (x2  +  y2  +  z2)1'2  [33] 

/z\ 

6  =  tan-^-^|  [34] 

y  cos  e  +  z  sin  e 

0  =  tan-  V - j  [35] 

16.  Due  to  the  wide  range  of  initial  particle  velocities  angles  y,  it 

was  decided  that  calculations  would  be  simplified  if  done  In  a  reference 
frame  aligned  at  y  to  Frame  1.  Consequently,  Frame  2  was  defined  for  each 
particle  test  such  that  the  starting  plane  was  perpendicular  to  the  initial 
particle  direction  (Figure  6).  All  transformed  into  this  frame  by  means  of 
the  following  rotation  matrix: 


II 

CM 

o 

cos 

Y 

0 

-sin  y 

0 

1 

0 

sin 

Y 

0 

COS  Y 

[36] 


17.  The  collection  envelope  for  this  sphere  geometry  was  necessarily 

two-dimensional.  Essentially,  the  half-interval  method  was  used  to 
calculate  z'  boundary  values  for  a  series  of  y  coordinates;  the  y‘ 
boundaries  themselves  were  also  calculated  by  means  of  the  half-interval 
method.  The  set  of  boundary  points  (y 1 ,  z‘)  were  then  connected  by  cubic 
splines  and  then  the  enclosed  envelope  area  calculated. 
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Target  Geometry  For  Horizontal  Cylinders 

18.  In  reference  Frame  1  (previously  defined  for  spherical  geometry) 

the  fluid  velocity  field  Is  given  by  the  following  equations5: 

17-  L2(X2  -  Z2)l 


U 


x 


=  U 


r4 


[37] 


Uy  =  0  [38] 

-2UxzL2 


where  r  =  (x2  +  z2)1/2  [40] 

It  is  readily  deduced  from  these  equations  that  the  particle  motion  will  be 
confined  to  the  x-z  plane  provided  that  the  Initial  particle  velocity  in 
the  y  direction  Is  zero.  That  velocity  was  set  to  zero  In  accordance  with 
the  previously  stated  initial  conditions;  therefore,  this  problem  was  two- 
dimensional  . 

19.  As  with  the  spherical  geometry.  Frame  2  was  defined  and  used  as 
the  main  frame  In  which  all  motion  and  envelope  calculations  were  made. 
Due  to  the  planar  motion,  the  envelope  was  only  one-dimensional,  requiring 
only  an  'upper'  and  a  'lower'  boundary  to  be  calculated  (Figure  7).  This 
case  Is  consequently  very  much  simpler  than  that  of  the  spheres. 

Target  Geometry  For  Vertical  Cylinders 

20.  The  theory  here  Is  essentially  the  same  as  in  the  horizontal 


UNCLASSIFIED 


UNCLASSIFIED 


/II 


cylinder  case.  The  only  change  Involves  the  flow  field  which  was  changed 
to  the  x-y  plane  from  the  x-z  plane  In  Frame  1  (Figure  8).  The  flow 
equations  are: 


U  =  U 
x 


1  - 


L2  (x2  -  y2) 


[40] 


-2UxyL2 


V  r1* 

(41] 

u2-o 

[42] 

where  r  =  (x2  +  y2) 1/2 

[43] 

21.  The  combination  of  an  Infinite  vertical  cylinder  with  the 

required  five  target  radii  starting  position  from  the  cylinder,  removed  the 
usefulness  of  Reference  Frame  2.  Therefore,  Frame  1  was  chosen  as  the 
frame  for  motion  and  envelope  calculations.  Due  to  the  motion  symmetry 

about  the  x-axIs  and  the  z-axis,  only  a  y-boundary  needed  to  be  calculated 

for  the  envelope;  clearly,  any  z-coordlnate  starting  position  will  result 
in  a  similar  trajectory  that  differs  only  by  a  fixed  z-directlon 
displacement. 

Target  Geometry  For  Man  Simulation 

22.  This  simulation  was  rather  crude,  and  it  comprised  many 

simplifying  approximations  (Figure  9).  First,  a  vertical  cylinder  of 
approximate  man  dimensions  was  used  since  the  determination  of  the  flow 
field  around  a  man  was  much  too  complicated  a  proposition  for  this  project. 
Second,  the  flow  field  of  an  infinitely  long  cylinder  was  used  because  a 
calculation  of  end  flow  conditions  around  a  finite  cylinder  was  deemed  too 
complicated  for  a  first  approximation  model  like  this  one.  Third,  the 
particles  which  impacted  on  the  top  of  the  cylinder  were  ignored. 
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23.  A  major  component  of  this  model  was  the  inclusion  of  a  wind 

gradient  to  mimic  the  Earth's  own  boundary  layer.  Essentially,  the 
horizontal  windspeed  is  a  function  of  height  above  the  ground.  The 
equation  used  was: 


[Z  \1'7 

U(Z)  =  Ul|-j  [44] 

The  base  point  (Zlt  Uj)  scales  the  curve;  for  this  application,  the  basis 
used  was: 


Z;  *  5.0  m 


Ux  =  1.0,  1.5,  2.5,  5.0  and  7.5  m/s 


[45] 

[46] 


Therefore,  five  tests  were  conducted  at  five  different  reference 
windspeeds. 


24.  The  flow  equations  for  vertical  cylinders,  equations  [40]  to  [42] 

were  modified  for  this  simulation  by  replacing  the  factor  U  with  as 
defined  in  equation  [44].  The  resulting  flow  equations  were: 


/  zy 

u*'u‘(  7J 


n 


1  - 


L2(x2  -  y2) 


[47] 


i 


[48] 


[49] 
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25.  As  for  the  vertical  cylinders,  Reference  Frame  1  is  used  for  all 
motion  and  envelope  calculations.  The  presence  of  the  wind  gradient 
destroys  the  problem  symmetry  in  the  z-direction;  therefore,  the  envelope 
on  the  starting  plane  was  necessarily  two-dimensional.  Its  calculation  was 
quite  similar  to  that  of  the  spherical  targets.  Specifically,  the  bounds 
in  the  z-direction  were  found  by  the  half  interval  method,  then  the  y- 
bounds  were  calculated  for  eleven  equally  spaced  z-coordinates  over  this 
Interval.  The  resulting  boundary  points  (z,y)  were  then  integrated  to 
yield  the  enclosed  envelope  area. 

26.  In  order  to  relate  the  chemical  concentration  on  a  man  to  that 
which  falls  on  the  ground,  an  extra  number  was  introduced,  termed  the 
ground  fraction  (6).  It  was  defined  by  equation  [50]: 

Man  Concentration  *  Ground  Concentration  x  G  [50] 

(This  man  concentration  is  based  on  frontal  cross-sectional  area,  not 
surface  area.) 

Calculation  of  G  combined  the  factors  of  collection  efficiency  with  the 
relative  areas  of  Impact  regarding  the  target  and  the  corresponding  ground; 
the  equation  was: 


Collection  Efficiency 

G - [51] 

tan  a 

In  this  equation,  a  is  a  representative  trajectory  angle  of  all  particles 
in  the  envelope;  note  that  the  trajectories  are  not  straight  lines  because 
of  the  wind  gradient.  The  angle  a  is  that  between  the  horizontal  and  the 
line  connecting  the  top  of  the  'man'  to  the  top  of  the  envelope. 
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Computer  Programs 

27.  During  the  course  of  the  project,  several  computer  programs  were 
written.  Four  collection  efficiency  programs  were  written;  one  for  each 
geometrical  case.  Two  programs  were  written  to  tabulate  and  plot  chemical 
particle  trajectories  for  the  sphere  and  horizontal  cylinder  cases;  these 
programs  were  used  primarily  for  verification  of  the  motion  model.  Final¬ 
ly,  a  number  of  utility  programs  were  written  to  tabulate  and  graph 
results.  All  programs  were  written  in  the  Honeywell  FORTRAN-77  language 
and  executed  on  the  Honeywell  CP-6  computer  at  DRES. 

28.  The  collection  efficiency  and  trajectory  plotting  programs  shared 
many  features,  most  Importantly,  the  initial  condition  and  motion 
calculations.  This  redundancy,  coupled  with  space  limitations,  has  limited 
the  program  listing  In  this  report  to  just  two  representative  programs: 
AEROSOL -8  which  performed  the  man  simulation  (Appendix  A),  and  AEROSOL-6, 
which  performed  the  trajectory  plotting  for  horizontal  cylinders  (Appendix 
B).  The  only  real  difference  between  these  programs  and  their  sister 
programs  was  in  geometry;  different  target  geometries  required  different 
envelope  calculations  and  sometimes  different  reference  frames,  as  has  been 
Illustrated  In  the  previous  section.  Implementation  of  the  background 
theory  on  the  computer  was  quite  similar  for  the  different  geometries,  and 
will  now  be  explained  in  detail  for  programs  AEROSOL-8  and  AEROSOL-6. 

29.  AEROSOL-8  was  designed  to  accommodate  five  sets  of  test 
conditions  (the  '700  loop',  commencing  line  50).  The  test  conditions, 
namely  wlndspeed,  target  radius  and  chemical  particle  density,  were 
obtained  from  the  data  file  AER0INF0.  The  program  then  proceeded  to 
calculate  collection  efficiencies  for  a  range  of  particle  sizes  which  were 
stored  in  the  array  SIZE(10).  This  loop  (the  '650  loop',  commencing  line 
63)  comprised  four  main  sections:  calculation  of  Initial  conditions  (lines 
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68-83),  calculation  of  the  z-dl recti  on  envelope  boundaries  (lines  86-138), 
calculation  of  the  y-dlrectlon  envelope  boundaries  (lines  141-161)  and 
calculation  of  resulting  collection  efficiency  (lines  164-172).  The 
Implementation  of  the  background  theory  In  these  sections  was  relatively 
stral ghtforward  except  for  the  following  points.  First,  In  order  to 
Implement  the  half-interval  method,  a  starting  position  which  results  in 
particle  Impact  must  be  found;  In  the  z-direction.  It  was  searched  for 
(the  *200  loop',  commencing  line  95),  but  In  the  y-direction  It  was  assumed 
that  y=0  resulted  in  impact  since  It  lay  on  the  target  centerline  and  would 
experience  no  sideways  drag  force.  Secondly,  preliminary  calculations 
showed  that  the  y-bounds  increased  monotonically  with  increasing  z- 
coordlnate,  therefore  the  previous  y-boundary  value  (variable  YMEM,  line 
143)  was  used  as  the  'Inside'  position  for  the  ensuing  half-interval 
calculation.  Finally,  note  that  the  problem  was  symmetrical  about  the  y- 
axis,  requiring  that  only  positive  y-boundarles  be  calculated  for  the 
envelope. 


30.  The  TRAJECTORY  subroutine  comprised  all  of  the  motion  calcula¬ 

tions  from  Initial  conditions  to  a  determination  of  particle  impact  or 
miss.  Upon  receipt  of  the  initial  conditions,  the  subroutine  sets  up  an 
iterative  loop  (the  '10  loop',  commencing  on  line  244)  to  'move'  the 
particle  towards  the  target.  The  Iteration  Involved  five  main  steps: 
calculation  of  local  flow  velocity  (lines  248-252),  calculation  of 
Reynold's  number  and  Drag  Coefficient  (lines  254-279),  solution  of  the 
differential  equations  of  motion  over  a  predetermined  time  Interval  (lines 
284-291),  testing  for  impact  or  miss  (lines  291-297)  and  adjustment  of  the 
differential  equation  step-size  If  the  Iteration  Is  to  continue  (lines  303- 
305).  Note  that  the  step-size  is  governed  by  the  choice  of  the  variable 
DT,  not  the  IMSL  routine  DVERK,  because  the  positional  dependence  of  the 
flow  velocity  requires  frequent  updating.  DVERK  could  not  accommodate  non¬ 
constant  coefficients,  and  in  practice  performed  only  one  step  per  call  due 
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to  the  choice  of  DT  which  determined  the  end  condition  TEND.  This  also 
meant  that  the  DVERK  accuracy  parameter  TOL  was  of  no  practical  Importance; 
Its  assigned  value  of  0.01  was  arbitrary. 

31.  The  remaining  two  program  subroutines  require  no  detailed 
explanation  above  the  internal  documentation  notes.  Note  that  a  partial 
output  is  Included  with  the  program  listing  in  Appendix  A. 

32.  AEROSOL-6  was  designed  to  plot  particle  'streamlines'  near  a 
horizontal  cylindrical  target.  It  was  restricted  to  one  size  of  particle, 
which  can  be  arbitrarily  set,  and  one  set  of  test  conditions.  The 
calculation  of  Initial  conditions  (lines  64-84)  was  identical  to  that  of 
AEROSOL-8,  with  the  exception  of  the  ability  to  perform  in  a  no  gravity 
environment;  hence,  this  program  could  recognize  and  plot  streamlines  for  a 
hypothetical  no  gravity  situation.  The  motion  analysis  section  (the  '60 
loop',  commencing  on  line  105)  was  the  same  as  that  of  AEROSOL-8. 

33.  The  one  aspect  of  AEROSOL-6  that  differs  from  AEROSOL-8  was  the 
trajectory  plotting.  To  fully  understand  the  details  of  this  plotting,  one 
must  study  the  CALCOMP  Electromechanical  Plotters  User's  Manual  (From 
California  Computer  Products,  Inc.).  As  far  as  this  project  was  concerned, 
CALCOMP  provided  a  list  of  subroutines  which  could  be  called  upon  to  draw 
graphs;  separate  subroutines  could  draw  axes,  plot  lines,  print  titles  and 
draw  other  graph  elements.  In  AEROSOL-6,  the  arrays  XVAL(500)  and 
ZVAL( 500)  were  used  to  store  the  x  and  z  coordinates  of  each  successive 
particle  location.  This  data  was  then  sent  to  the  CALCOMP  System  and 
plotted  (lines  211-241).  Examples  of  these  trajectory  plots  as  given  in 
Figures  10  and  11,  and  a  sample  output  is  included  in  Appendix  B  with  the 
program  listing. 
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RESULTS 


34.  The  computer  generated  collection  efficiency  results  for  all  four 

geometrical  cases  are  listed  In  Tables  I-VIII.  Figures  10  and  11  show 
trajectory  pictures  generated  by  AEROSOL-6.  Figures  12-15  show 
representative  and  comparative  graphs  based  on  the  data  contained  In  Tables 
I-VIII.  Although  a  detailed  analysis  of  these  results  will  be  conducted  In 

the  next  section,  the  acquisition  and  content  of  these  results  require  some 

explanations.  Note  that  all  results  are  In  MKS  units  unless  otherwise 
stated. 

35.  Sphere  runs  SI  to  S10,  and  horizontal  cylinder  results  Cl  to  CIO 

were  all  performed  under  gravity  and  no  gravity  conditions  to  calculate  the 
effect  of  gravity  on  the  problem.  These  results  are  listed  side  by  side  In 

Tables  I,  II,  IY  and  V.  At  the  bottom  of  each  double  column  are  two  sets 

of  numbers  labelled  'maximum  positive  change',  and  'maximum  negative 
change'.  These  are  simply  the  greatest  divergence  values  between  the  two 
sets  of  data,  recorded  as  either  a  positive  or  a  negative  divergence 
relative  to  the  gravity  values.  The  numbers  In  brackets  are  the  particle 
sizes  corresponding  to  those  divergences. 

36.  Because  spheres  were  the  first  geometrical  case  to  be  studied,  a 
few  extra  runs  were  conducted  In  order  to  evaluate  the  validity  and  the 
accuracy  of  the  motion  model.  The  accuracy  test,  run  SI  1 ,  used  half  the 
step-size  as  run  S2  In  order  to  check  the  numerical  accuracy  of  the  method. 
In  actuality,  several  accuracy  checks  were  conducted,  resulting  in  many 
step-size  modifications  until  the  final  accuracy  level  was  achieved.  The 
similarity  test,  runs  S12  and  S13,  was  conducted  in  an  attempt  to  validate 
the  motion  model.  According  to  the  principles  of  fluid  mechanics,  tests 
with  the  same  set  of  non-dimensional  numbers  must  yield  identical  results. 
Here,  there  are  three  non-dimensional  numbers  which  describe  the  problem: 
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CQ>  ReQ,  F.  These  were  kept  constant  for  runs  S2,  S12,  and  S13,  but  the 
constitutive  parameters  {U,  L,  F,  P,  g)  were  varied  to  learn  whether  or  not 
the  results  would  change.  A  change  In  the  results  would  Indicate  that  the 
motion  model  was  flawed. 

37.  The  computer  time  required  for  solution  of  test  runs  varied 
considerably  with  the  problem  geometry.  The  one-dimensional  nature  of  the 
horizontal  and  vertical  cylinder  envelopes  resulted  in  very  little  computer 
execution  time,  typically  ten  minutes  or  less  per  test.  (By  way  of 
clarification,  a  'test'  refers  to  the  calculation  of  collection 
efficiencies  for  all  ten  particle  sizes  under  a  given  set  of  windspeed, 
target  size  and  particle  density  values.  Each  column  in  Tables  I-VIII 
represents  one  test.)  The  two-dimensional  envelopes  for  the  sphere  cases 
typically  required  one  and  a  half  hours  per  test.  The  man  simulation  was 
somewhat  peculiar.  Because  the  particles  travelled  almost  horizontally  in 
test  Ml,  it  required  only  twenty  minutes;  but  the  more  vertical 
trajectories  of  test  M4  and  M5  resulted  in  much  longer  trajectories  (note 
that  the  horizontal  distance  travelled  remained  constant  so  as  to  minimize 
the  cylinder  flow  perturbation  at  start)  and  execution  times  of  up  to  six 
hours  per  test.  Actually,  the  collection  efficiencies  for  the  2000  un 
particle  were  not  calculated  for  tests  M4  and  M5,  and  the  1000  un  particle 
collection  efficiency  was  not  calculated  for  M5.  The  reason  was  that  hours 
of  execution  time  would  have  been  required  for  each  particle,  a  cost  which 
was  not  thought  to  be  worthwhile.  The  absence  of  these  three  values  is 
indicated  in  Tables  VII  and  VIII  by  a  negative  value. 

38.  Finally,  note  that  the  graphs  plotted  in  Figures  12-15  were  done 
using  small  plotting  programs  written  during  the  project.  These  programs 
are  not  listed  in  this  report. 
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DISCUSSION 


39.  The  computed  similarity  and  accuracy  test  results  (Table  III) 
verified  the  motion  model  and  provided  an  Indication  of  its  accuracy.  The 
accuracy  test  clearly  showed  that  numerical  accuracy  Improved  with 
Increasing  particle  size.  This  was  because  larger  particles  were  less 
affected  by  drag  forces,  and  it  was  the  drag  force  which  exhibited  non¬ 
linear  behaviour  with  position,  making  It  the  most  difficult  aspect  to 
numerically  approximate.  This  also  accounted  for  the  observation  that  all 
errors  were  positive,  because  the  total  Influence  of  the  drag  force  (and 
hence  the  greatest  particle  deflection)  will  only  be  attained  in  the  limit 
as  the  number  of  steps  approaches  Infinity.  Hence,  the  drag  force  was 
underestimated  by  this  model.  Since  the  overall  collection  efficiency  was 
the  chief  result  sought,  the  relative  error  was  of  less  Importance  here; 
that  is  to  say,  the  large  relative  error  of  the  smallest  particles  was  made 
Insignificant  by  the  very  small  collection  efficiencies  Involved.  Based  on 
the  absolute  errors  tabulated,  the  maximum  error  present  In  the 
calculations  was  on  the  order  of  +  0.003  for  the  collection  efficiency  (+ 
0.3%  for  the  values  In  Tables  I-VII).  It  should  be  noted  that  the  computer 
model  worked  to  four  significant  digits  In  distance  values  (0.1  mm  or  100 
ym).  Since  the  targets  were  generally  0.1  m  in  radius,  0.1  n«n  represents 
an  accuracy  to  0.1%. 

40.  The  results  of  the  similarity  test  were  within  numerical  error 
for  all  three  runs,  provided  that  one  qualification  to  the  collection 
efficiency  error  value  of  +  0.003  be  accepted;  specifically,  that  the 
absolute  collection  error  was  a  function  of  the  target  size  In  addition  to 
step-size.  It  seems  plausible  that  the  relative  magnitude  of  the  step  size 
to  the  target  size  would  Influence  numerical  accuracy,  since  an  increase  In 
the  target  size  would  attenuate  the  rate  of  change  of  the  drag  force  over 
distance,  allowing  a  greater  accuracy  for  the  same  step  size.  As  proof. 
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note  that  the  S12  and  S13  results  are  within  .003  of  each  other,  with  the 
S13  values  uniformly  lower;  this  suggests  that  its  larger  target  radius 
( L=0 . 15 )  improved  numerical  accuracy,  since  the  step  sizes  were  equal. 
Conversely,  the  small  target  size  of  test  S12  seemed  to  have  degraded 
numerical  accuracy;  the  S12  values  were  up  to  .008  higher  than  the  S2 
values.  In  summary,  the  results  of  tests  S2,  S12  and  S13  were  deemed 
sufficiently  close  as  to  be  judged  the  same  to  within  numerical  error. 

41.  Figures  10  and  11  qualitatively  demonstrate  many  of  the  aspects 
of  this  motion  problem.  The  smaller  particle  (Figure  10)  showed  a  mostly 
horizontal  trajectory  which  indicated  low  terminal  velocity.  Close  to  the 
cylinder,  all  trajectories  were  substantially  deflected  to  the  extent  that 
only  one  particle  Impacted  on  the  cylinder.  This  was  an  Indication  that 
the  collection  efficiency  would  be  low  for  this  small  (50  jn)  particle. 
Figure  11  demonstrated  the  aspects  of  large  particle  motion.  The  trajec¬ 
tory  was  much  steeper  due  to  a  higher  terminal  velocity.  There  was  almost 
no  particle  deflection;  consequently,  one  would  expect  a  high  collection 
efficiency  for  this  large  (250  ym)  particle.  Reference  to  run  C4  which  had 
the  same  test  conditions  as  In  Figures  10  and  11,  yielded  the  expected 
magnitude  of  collection  efficiencies:  a  low  14.11%  for  50  un,  and  a  high 
93.70%  for  250  urn.  The  qualitative  model  verification  by  these  and  many 
other  trajectory  pictures  supported  the  similarity  tests  and  led  us  to 
conclude  that  the  motion  model  was  valid. 

42.  The  collection  efficiency  test  results  for  spheres  and  cylinders 
yielded  many  noteworthy  features,  most  of  which  will  be  discussed  in  the 
next  few  pages.  We  will  start  with  the  effect  of  gravity  on  the  motion 
results. 


43.  All  of  the  sphere  and  horizontal  cylinder  tests  (SI  -  S10  and  Cl 

-CIO)  produced  results  similar  to  that  shown  In  Figure  12.  The  sigmoidal 
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shape  of  the  curve  agreed  with  the  previous  theoretical  work.7  The  most 
notable  feature  of  the  double  curve  plot  in  Figure  12  was  the  crossing  of 
the  two  curves  between  100  pm  and  250  na  particle  sizes;  this  was  thought 
to  be  due  to  the  following  reasons.  For  larger  particles,  gravity  caused 
the  particles  to  fall  with  a  terminal  velocity  close  to  or  greater  than  the 
horizontal  free  stream  velocity.  Hence,  the  particle's  inertia  was 
significantly  greater,  rendering  the  particle  much  less  susceptible  to 
deflection  by  the  diverging  fluid  streamlines  around  the  target.  This 
resulted  in  a  collection  efficiency  with  gravity,  as  evidenced  by  the 
larger  particles  in  Figure  12. 

44.  The  gravity-decrease  effect  on  the  smaller  particles  was  much 

more  difficult  to  explain.  First,  the  terminal  velocity  was  low,  virtually 
insignificant  compared  to  the  free  stream  velocity;  therefore,  the 
particle's  Inertia  was  not  noticeably  higher.  This  allowed  a  second  factor 
to  make  a  discernible  Impact  on  the  motion;  this  factor  was  a  motion 

asymmetry  around  the  target  due  to  gravitational  effects.  The  envelope 
results  for  test  C3  (Table  9)  Illustrate  this  asymmetry.  It  can  be  seen 
that  gravitational  effects  decrease  the  upper  bound  more  than  they  Increase 
the  lower  bound  for  particles  of  50  -  100  pm,  which  was  the  range  of  the 
gravity  decrease  effect  In  Figure  12.  The  explanation  of  this  asymmetry 
was  as  follows.  Near  the  upper  boundary,  the  fluid  streamlines  are 
deflected  upwards  by  the  target,  in  effect  flowing  crossways  to  the 
particle  motion  vector,  thereby  increasing  the  local  Reynold's  numbers  and 
hence  the  drag  force,  resulting  in  a  smaller  boundary.  At  the  lower 

boundary,  however,  the  fluid  streamlines  are  deflected  in  the  direction  of 
the  particle  motion  vector  (that  is,  diagonally  downwards)  thereby  reducing 
the  local  Reynold's  number  and  the  drag  force,  resulting  in  a  larger 

boundary.  The  non-linear  nature  of  the  problem  was  such  that  the  former 

effect  was  greater  than  the  latter  effect,  resulting  in  lower  collection 
efficiencies  for  particles  in  this  motion  regime  when  gravity  was  Included. 


UNCLASSIFIED 


UNCLASSIFIED 


/22 


45.  Generally,  gravitational  effects  altered  collection  efficiencies 
In  these  tests  by  less  than  ±3%  (Tables  I,  II,  IV,  V).  The  exceptions  were 
tests  Involving  low  wlndspeeds,  speeds  of  1.5  m/s  or  less.  The  greatest 
difference  was  found  In  test  S5;  there,  gravity  added  21.2%  to  the 
collection  efficiency  of  100  pm  particles.  Note  that  the  crossover  point 
between  the  gravity  and  no  gravity  curves  In  this  test  was  not  present;  all 
differences  were  positive.  Generally,  the  crossover  point  decreased  in 
particle  size  with  decreased  windspeed. 

46.  In  summary,  the  gravity  effect  on  collection  efficiencies  was 
negligible  except  in  cases  of  low  windspeed.  Other  facets  of  the  sphere 
and  cylinder  studies  will  now  be  explored. 

47.  Figure  13  compared  the  collection  efficiencies  of  the  vertical 
cylinder,  the  horizontal  cylinder  and  cylinder  no  gravity  cases,  under  the 
same  test  conditions.  (Note  that  in  the  absence  of  gravity,  the  vertical 
and  horizontal  cylinders  were  geometrically  equivalent.)  The  crossover  of 
the  no  gravity  and  horizontal  cylinder  curves  occurred  at  500  tin;  the 
gravity  effect  was  minimal  here  as  would  be  expected  from  the  high 
windspeed.  The  vertical  cylinder  curve  was  lower  than  either  of  the  other 
curves,  except  for  a  brief  particle  range  around  100  urn.  The  reduced 
collection  efficiency  was  easily  explained:  as  particles  approached  the 
cylinder,  they  initially  deccelerated,  thereby  increasing  the  trajectory 
angle  relative  to  the  horizontal  (the  terminal  velocity  remains  almost 
constant)  and  provided  more  time  for  the  diverging  fluid  streamlines  to 
deflect  the  particle.  The  exception  at  100  urn  was  probably  a  result  of  the 
asymmetrical  particle  flow  around  the  horizontal  cylinder,  postulated 
before  as  the  explanation  for  the  gravity-decrease  effect.  The  vertical 
cylinder  possessed  symmetrical  flow  conditions  and  was  therefore  not 
affected.  However,  the  steepening  trajectory  effect  with  vertical 
cylinders  still  lowered  the  collection  efficiency  relative  to  the  no 
gravity  case;  the  reduction  was  just  marginally  less  than  that  of  the 
horizontal  cylinders. 
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48.  Figure  14  compared  the  two  cylinder  geometries  to  spheres,  under 

slightly  different  test  conditions  than  In  Figure  13.  Spheres  clearly 
possessed  higher  collection  efficiencies  than  cylinders  for  all  particle 
sizes.  This  can  be  understood  in  light  of  the  fact  that  cylinders  perturb 
the  flow  more  than  spheres.  In  the  sense  that  cylinders  represented  a 
greater  obstacle  to  the  flow  and  thus  cause  faster  fluid  motion  around  the 
periphery.  This  larger  fluid  velocity  represented  a  greater  drag  force 
which  tended  to  deflect  the  particles  away  from  the  target;  hence,  the 

collection  efficiencies  were  lower.  The  relationship  between  the  vertical 
and  horizontal  cylinder  curves  was  the  same  as  in  Figure  13  except  that  the 
difference  here  was  smaller  because  of  the  smaller  target  size.  It  should 
be  noted  that  Figures  13  and  14  were  representati ve  of  all  of  the  test 
results  listed  in  Tables  I  to  VI,  and  were  not  just  the  product  of  those 
specific  test  conditions. 

49.  Some  general  observations  on  the  sphere  and  cylinder  studies  will 

now  be  made.  In  all  three  cases,  a  decrease  in  windspeed  resulted  in 
reduced  collection  efficiency  for  any  given  particle  size.  Evidently,  the 
reduced  drag  force  was  more  than  compensated  for  by  the  decrease  in 

particle  inertia.  Increased  particle  density  appeared  to  merely  shift 
collection  efficiency  values  from  larger  to  smaller  particles; 

equivalently,  plots  of  collection  efficiency  vs  log  (particle  radius)  were 
translated  left.  The  reason  for  this  was  that  particle  terminal  velocity 
and  inertia  were  increased,  rendering  the  particles  more  difficult  to 
deflect.  A  reduction  in  target  size  had  the  same  effect  as  increased 
particle  density.  The  explanation  in  this  case,  however,  was  that  larger 
targets  create  more  far-reaching  flow  perturbations,  the  net  effect  of 
which  was  subject  to  incoming  particles  to  greater  deflecting  drag  forces; 
conversely,  smaller  targets  resulted  in  less  deflecting  drag  forces. 

50.  The  remaining  discussion  will  focus  on  the  man  simulation 

results.  Although  the  ground  fraction  values  are  of  most  importance  here, 
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a  brief  comparison  of  the  collection  efficiency  results  In  Tables  VI  and 
VII  needs  to  be  done.  The  only  difference  between  the  two  tests  was  that 
the  'M'  tests  (Table  VII)  incorporated  a  velocity  gradient  from  the  ground 
up.  This  velocity  gradient  significantly  lowered  the  collection  efficiency 
for  all  particle  sizes.  The  reason  was  that  as  the  particles  fell,  the 
wlndspeed  decreased  which  In  turn  decreased  the  particle's  horizontal 
speed.  As  In  all  previous  cases  of  reduced  wlndspeed,  this  must  result  In 
reduced  collection  efficiency. 

51.  The  ground  fraction  results  demonstrated  several  noteworthy 

features.  Most  striking  were  the  greater  than  unity  ground  fractions  for 
some  of  the  particles  In  high  wlndspeed  tests.  This  was  due  to  almost 
horizontal  particle  trajectories  for  these  cases;  near  horizontal 
trajectories  will  result  In  low  ground  concentration  since  the  particles 
will  be  distributed  over  a  large  area.  Clearly,  objects  standing 

vertically  In  such  a  situation  could  receive  greater  concentrations  than 
the  ground. 

52.  Figure  15  showed  three  curves  corresponding  to  tests  Ml,  M2  and 

M3.  The  central  peak  In  each  curve  resulted  from  the  Interaction  of  two 
effects.  For  small  particles,  the  collection  efficiency  was  so  small  that 
the  ground  fraction  was  zero;  for  large  particles,  their  near  vertical 
trajectories  meant  that  relatively  few  could  impact  on  the  vertical  sides 
of  the  cylinder,  so  that  the  ground  fraction  was  again  near  zero.  Both 
ground  fraction  reducing  effects  decreased  toward  the  opposite  end  of 
the  particle  size  spectrum;  therefore,  the  largest  ground  fractions 

occurred  in  the  middle  region,  in  which  neither  effect  dominated.  This 
peak  migrated  from  50  pm  in  the  upper  curve  to  100  nn  In  the  lower  curve. 

53.  The  tremendous  effect  of  wlndspeed  was  well  Illustrated  by  Figure 

15.  The  lower  wlndspeeds  possessed  very  low  ground  fractions;  note  that  a 
maximum  ground  fraction  of  0.17  was  calculated  for  the  1.0  m/s  wlndspeed 
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case.  The  reason  for  this  was  that  particle  trajectories  were  near 
vertical  for  such  low  windspeeds;  hence,  little  particle  impaction  could 
occur  on  the  vertical  cylinder  sides. 

54.  The  man  simulation  model  in  this  project  was  undoubtedly  crude. 
Nevertheless,  the  essential  aspects  of  this  kind  of  man  simulation  problem 
were  believed  to  have  been  demonstrated,  even  though  the  flow  geometry  was 
drastically  simplified.  Some  speculative  conclusions  will  now  be  drawn 
from  the  ground  fraction  data. 

55.  The  curves  of  Figure  15  indicate  that  the  largest  aerosol 
particles  are  not  suitable  for  impacting  on  a  standing  man;  in  fact,  the 
best  particle  size  is  around  100  pm.  This  must  be  viewed  in  light  of  two 
important  qualifiers:  first,  the  smaller  particles  will  travel  further  from 
the  dissemination  point,  resulting  in  lower  ground  concentrations  to  begin 
with;  and  second,  aerosols  from  materials  with  some  volatility  will  tend  to 
evaporate  as  they  move  downwind,  so  that  the  smaller  particles  might 
disappear  altogether.  Note  that  the  evaporative  characteristics  will  also 
influence  ground  persistence  of  the  chemical,  which  is  another  vital 
consideration.  Nonetheless,  these  ground  fraction  results  suggest  that  an 
upper  limit  for  ideal  particle  size  for  impaction  on  a  man  may  exist. 

56.  The  subject  of  man  motion  under  this  kind  of  chemical  particle 
bombardment  was  not  considered  in  the  project.  Although  authoritative 
comments  will  have  to  wait  until  detailed  work  is  done,  there  is  one 
speculation  that  needs  to  be  recorded  here.  Specifically,  if  the  man  were 
walking  in  the  direction  of  the  wind,  the  particle  trajectories  would 
assume  a  more  vertical  shape  in  the  man's  frame  of  reference.  This  would 
be  equivalent  to  a  man-stationary,  reduced  windspeed  problem  such  as  was 
studied  in  this  project;  and  according  to  those  results,  the  ground 
fraction,  and  hence  the  man  contamination,  would  be  reduced.  Incorporation 
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of  a  moving  target  Into  this  simulation  model  Is  a  logical  next  step  for 
research,  one  that  would  help  to  resolve  the  speculation  suggested  above. 

CONCLUSIONS 


57.  The  numerical  tests  conducted  with  spherical  cylindrical  targets 
Indicated  that  gravitational  effects  altered  the  collection  efficiencies 
insignificantly,  on  the  order  of  ±  3%  (absolute),  provided  that  the  free 
stream  velocity  was  2.5  m/s  or  higher.  Lower  free  stream  velocities 
resulted  In  much  greater  gravitational  effects,  up  to  21.2%  (absolute)  for 
100  pm  particles  Impacting  on  a  spherical  target.  Changes  In  the  particle 
density  and  target  size  were  found  to  have  a  negligible  effect  on  the 
Importance  of  gravity  In  the  problem. 

58.  The  man  simulation  tests  Indicated  that  a  man  would  receive  the 
most  chemical  relative  to  the  ground  concentration  for  particles  on  the 
order  of  100  pm  radius.  In  fact,  he  could  receive  up  to  6.9  times  the 
ground  concentration.  Although  there  were  mitigating  factors,  the  analysis 
suggested  that  an  upper  limit  may  exist  for  the  ideal  particle  size  In 
considering  Impaction  on  a  standing  man. 
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TABLE  III 

MISCELLANEOUS  SPHERE  RUNS 

» 

ACCURACY  TEST 


S2 

Sll 

5.0 

5.0 

0.1 

0.1 

25.5 

25.5 

1000 

1000 

9.81 

9.81 

A8S.  RE 
ERROR  ER 


PARTICLE 

RAOIUS 


.00044 

.00640 

.11263 

.44986 

.63808 

.74054 

.92237 

.97778 

.99656 

1.00032 


.00032 

.00515 

.11008 

.44751 

.63635 

.73904 

.92182 

.97771 

.99656 

1.00032 


.00173 

.00150 

.00060 

.00007 


37.5% 

24.3% 

2.3% 

0.5% 

0.3% 

0.2% 

0.06% 

0.007% 


SIMILARITY  TEST 


S2 
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5.0 

5.0 

5.0 

0.1 

0.05 

0.15 
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25.5 

25.5 

1000 
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1500 

9.81 

19.62 

6.54 

PARTICLE 

RADIUS 
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TABLE  IX 

ENVELOPE  RESULTS  FOR  TEST  C3 
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RADIUS 

(microns) 
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BOUND 

(m) 

LOWER 

BOUND 

(m) 

NO  GRAVITY 

BOUND 

(m) 

10 

±.0001 

15 

.0001 

25 

-.0006 

.0016 

50 

+.0162 

-.0299 

.0236 

75 

+.0314 

-.0523 

.0439 

100 

+.0435 

-.0668 

.0579 

250 

+.0827 

-.0967 

.0873 

500 

+  .0971 

.0959 

1000 

+  .1004 

-.1009 

.0994 

2000 

+  .1019 

-.1019 

.1014 

The  boundary  values  listed  here  are  z'  co-ordinates  (Frame  2).  The  target 
was  0.1  m  in  radius. 
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Figure  5 

SPHERE  REFERENCE  FRAMES  0  AND  1 
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Figure  8 

VERTICAL  CYLINDER  GEOMETRY 
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Figure  9 

MAN  SIMULATION  GEOMETRY 
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FIGURE  10 

PARTICLE  S  l  ZE  =  50-  MICRONS  W|NDSPEED=  1 . 5M/S 
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APPENDIX  A 

LISTING  OF  PROGRAM  AEROSOL-8  (MAN  SIMULATION) 
SAMPLE  OUTPUT  FOR 
FIRST  SET  OF  TEST  CONDITIONS  (RUN  1) 
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